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ABSTRACT 

We present the detection and characterization of the transiting warm Jupiter KOI-12b, first identified with Kepler 
with an orbital period of 17.86 days. We combine the analysis of Kepler photometry with Doppler spectroscopy and 
line-profile tomography of time-series spectra obtained with the SOPHIE spectrograph to establish its planetary nature 
and derive its properties. To derive reliable estimates for the uncertainties on the tomographic model parameters, 
we devised an empirical method to calculate statistically independent error bars on the time-series spectra. KOI-12b 
has a radius of 1.43±0.13 itj up and a 3<r upper mass limit of 10Mj up . It orbits a fast-rotating star (usini* = 
60.0±0.9kms -1 ) with mass and radius of 1.45±0.09Mg and 1.63±0.15 Rq, located at 426±40pc from the Earth. 
Doppler tomography allowed a higher precision on the obliquity to be reached by comparison with the analysis of 
the Rossiter-McLaughlin radial velocity anomaly, and we found that KOI-12b lies on a prograde, slightly misaligned 

_|_3 Q° 

orbit with a low sky-projected obliquity A = 12 . 6 - 2)9 . The properties of this planetary system, with a 11.4 magnitude 
host-star, make of KOI-12b a precious target for future atmospheric characterization. 


Key words, planetary systems - Stars: individual: KOI-12 


1. Introduction 

The transit of an exoplanet across its rotating host star 
distorts the apparent stellar line shape by removing the 
part of the profile emitted by the occulted portion of the 
star. This induces anomalous variations in the measured 
stellar radial velocity during the tr ansit, know n as the 
Rossiter-McLaughlin (RM) anomaly (Holt 1893 Rossiter 


1 1924 McLaughlin 1924). The shape of the anomaly as a 


function of time depends on the value of the sky-projected 
spin-orbit obliquity A, which is the angle in the plane of the 
sky between the projection of the stellar spin axis and the 
projection of the orbital angular momentum vector. Most 
obliquity measurements have been obtained using Doppler 
spectroscopy, but complementary techniques also make use 


et al. 2011 

Sanchis-Ojeda et al.| 2012), gravity darkening 

(Szabo et a 

. 2011 Barnes et al. 2013 

), asteroseismology 
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( Chaplin et al. 2013 ) and Doppler tomography Collier 
Cameron et al.||2010al. 


Obliquity is a good tracer of planetary system histories, as 
different formation scenarios result in different spin-orbit 
angles. About thirty misaligned systems (|A| >30° and in¬ 
consistent with A =0°) have bee n identified today over more 
than eighty measured system^ 4 ] (|Albrecht et al.|2012 Crida 
fe Batygin|2014), including some with retrograde or nearly 


polar orbits (e.g'|Winn et aL|2009[|Narita et al. 2010 Triaud 

et al. |2010 Hebrard et al. |2011|). Most of these measure¬ 
ments have been done for cfose-in, isolated giant planets. 
While it is commonly accepted that hot-Jupiters form be¬ 
yond the snow line and later migrate toward the star, many 
unknowns remain about how the migration occurs. Their 
wide distribution of obliquities favours misaligning scenar¬ 
ios where massive giant planets have been brought in by 
planet-planet (or planet-star) scattering, and Kozai migra- 


1 the Holt-Rossiter-McLaughlin Encyclopaedia: 
http://www.physics.mcmaster.ca/~rheller/ 
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tion with tidal friction 

(see e.g. Fabrycky & Tremainej2007 

Guillochon et al. 2011 

Naoz et al. 2011). Some models show 


tained through its interactions with the disk (Teyssandier 


et al. 2013). Instead, more standard scenarios implying disk 


migration are expected to conserve the initial alignment 
between the angular mo men t ums of the disk and of the 
planetary orbits (see, e.g., Lin et al. JT996). Alternatively, it 
has been proposed that the orbit still reflects the orienta¬ 
tion of the disk, with the stellar spin instead having moved 
away, eithe r early-on t hrough magnetosphere-disk interac¬ 
tions (Lai et ah 2011) or later through elliptical tidal in¬ 
stability ( Cebron et al.|20lT ). To understand whether such 
scenarios and the resulting wide distribution of obliquities 
are specific to massive close-in exoplanets, it is necessary 
to sample all types of planetary systems. While obliquity 
measurements have now extended to smaller planets (e.g. 


Winn et al.|2010b Hirano et al.||2012 Albrecht et al.|2013| 

Bourrier fe Hebrard||2014 Lopez-Morales et al.|2014 ), the 


ev 


remain little known for long-period exoplanets (live case^] 
with P ^11 days) or those orbiting fast-rotating stars (five 
case+Mwith wsini* ^20kms _1 ). 

The Kepler candidate KOI-12.01 (Table |T|) offers the op¬ 
portunity to probe spin-orbit angles in both domains, as 



the planetary nature of the transit signal based on Kepler 
photometry alone. Binary systems can mimic a planetary 
transit signal and are the so urces for a non negligible part 
of the Kep ler candidates (e.g. Santerne et al.|2012a Fressin 


et al.|2013 ). This is particularly true for candidates showing 
transit depths around 1% and thus possibly corre spo nding 
to giant plan ets, which is the case for KOI-12.01. Santerne 


et al. (2012a) further investigated the nature of KOI-12.01 
using radial velocity follow-up, but it proved unsufficient to 
solve the candidate. In recent years, Doppler tomography 
has been used to study the alignement and properties of/a 


growing number of planetary systems (HD 189733, Collier 
ICameron et al.|2010a| WA SP-3, Miller et al.|2010| CoRoT- 


11b, 


2012 


Gandolfi et al.||20l2t WA SP-32, -38, -40, Brown et al.| 


Kelt-7b, Bieryla et al. 2015) but also to assess the 


planetary nature of the transiting gas giant WASP-33b 


(Collier Cameron et al. 

2010b). This technique, similar to 

the methodology initial 

y developped to model the RM ef- 

feet in binary stars ( 

Albrecht et al. 

2007 

Albrecht et al. 


into its different components, namely the stellar and instru¬ 
mental profile, the limb-darkened rotation profile, and the 
travelling signature caused by the transiting planet. The 
large rotational broadening of the stellar lines for stars like 
WASP-33 and KOI-12 introduces errors when fitting the 
radial v elocities of the RM anomaly with anal ytical formu¬ 
lae (e.g. Triaud et al.|2009 Hirano et al.|2010 ). In contrast, 


Kepler-30 d, Kepler-30 c, Kepler-30 b, Sanchis-O.jeda et al. 


d2012|) ; HD 80606 b, |Hebrard et al. | (|2010|); HD 17156 b, |Narita 
et af.| (|2008 1 
v 


HAT-P-32 b and HAT-P-2 b. 


Albrecht et al. 


(2010b|; Keit-7b, Bieryla et al. (20151 


( 2012 ); 


CoRoT- 


11 b, | Gandolfi et al.|(fMT2|: WASP-33 b JColller Cameron et al. 


Doppler tomography is best-suited to the analysis of such 
systems, as it models directly the missing starlight signa¬ 
ture caused by the planet and reponsible for the radial ve¬ 
locity anomaly. 

Here we combine Kepler photometry and SOPHIE time- 
series spectroscopy to assess the nature of KOI-12.01 and 
determine the planetary system properties. Observations 
and data reduction are described in Sect. [2j The analysis 
of Kepler photometry is presented in Sect. [4j Radial ve¬ 
locities are analyzed in Sect. [5j and spectroscopic data are 
also used in Sect [6] to perform line-profile tomography. The 
planetary nature of KOI-12.01 is discussed in Sect. [TJ and 
we summarize our results in Sect. |U 


Table 1 . IDs, coordinates and magnitudes of KOI-12. 


Parameter 

Value 

Kepler object of interest 

KOI-12 

Kepler exoplanet catalog 

Kepler-? 

Kepler Input catalog 

KIC 5812701 

2MASS ID 

19494889+4100395 

WISE ID 

194948.89+410039.6 

RA (2000.0) 

19' 1 49" m 48 s .90 

Dec. (2000.0) 

+41°0 39".56 

Kepler magnitude 

11.353 

SDSS-G 

11.571 

SDSS-R 

11.280 

SDSS-I 

11.245 

2MASS-J 

10.461+0.020 

2MASS-H 

10.287+0.022 

2MASS-K, 

10.234+0.018 

WISE-W1 

10.189+0.023 

WISE-W2 

10.198+0.020 

WISE-W3 

10.015+0.042 


2. Observations and data reduction 


2.1. Photometric data with Kepler 
Seventy-four transits of the planetary candidate KOI-12.01 


were observed with Kepler (Batalha et al. 2011 2013), 


with a period of 17.9 days and transit depth of about 1 %. 
No transits with different periods were detected in any of 
the light curves, and there are thus no signs that KOI- 
12 is a multiple transiting system. The Kepler photometry 
was acquired in long-cadence (LC) and short-cadence (SC) 
data, with one point per 29.42 minutes and 58 seconds, 
respectively. We used in the present analysis a combina¬ 
tion of data in both cadences, with 28 LC transits and 46 
SC transits (Table [2]) . We used the light-curve of quarters 
Q0 to Q17 reducea by the Photometric Aperture Kepler 
pipeline, that accounts for barycentric, cosmic ray , back- 
ground, and so-called argabrightening corrections (Jenkins 
et al.| 2010), publicly available from the MAST archive 
(http://archive .stsci.edu/kepler). 

In addition to the transits, the Kepler light curve shows 
quasi-periodic photometric variability with typical ampli¬ 
tudes of 1000 ppm on timescales similar to the transit du¬ 
ration (~7.5h). They are probably due to inhomogeneities 
of the stellar surface (spots, plages, etc) modulated with 
the rotation of the star. Before modeling the transits, we 
normalized fragments of the light curves by fitting an out- 
of-transit 2 nd order polynomial. Because of the relatively 
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short timescale of the stellar variability, we limited this fit 
to ~1.9 hours of out-of-transit data (corresponding to more 
than 200 photometric measurements in the SC data and 
about 8 in the LC data). Given the high quality of the 
Kepler LC photometry, we found that this was sufficient to 
detrend all but one of the transits. While we removed most 
of the stellar variability with a timescale greater than the 
transit duration, stellar variability with a shorter timescale 
is not corrected and may explain the small increase in the 
scatter of the photometric data observed during the tran¬ 
sit in Fig. [l] Planet-spot occultations might also explain 
this increased scatter, which does not affect significantly 
the derived parameters but may slightly increase their un¬ 
certainties. Since the Kepler spacecraft rotates four times 
a year, the crowding values are different between seasons. 
We thus produced four crowding-uncorrected de-trended 
light curves for each season, using the crowding estimates 
provided by the MAST database (given the relative bright¬ 
ness of the target, we estimated an uncertainty of 0.2%). 
This allowed us to account for differential crowding values, 
noises, and out-of-transit fluxes in the transit modeling. 
Fig- [D shows the corresponding phase-folded transit light 
curves, once normalized and corrected for crowding. 



Orbital phase 

Fig. 1. Long-cadence transit light curve of KOI-12.01, nor¬ 
malized and phase-folded over the Kepler observations. The 
best fit to the data is displayed as a red line, with resid¬ 
uals in the lower panel. One of the light curves was found 
to be abnormally deep, and there is also a small increase 
in the residual scatter during the transits. These variations 
should not affect the derived parameters, but may slightly 
increase their uncertainties. 


2.2. Spectroscopic data with SOPHIE 


Spectroscopic measurements of KOI-12 were made with the 
spectrograph SOPHIE at the 1.93-m telescope of Haute- 
Provence Observatory (France). Two transits of KOI-12.01 
were observed in June 2012 (27 measurements, 20 dur¬ 
ing the transit) and June 2013 (18 measurements, 12 dur¬ 
ing the transit), in order to detect the RM anomaly and 
perform line-profile tomography. Additional measurements 
were taken far from the transit in 2011 (2 in February and 
June 2011; Santerne et al. 2012a) and in 2014 (2 in May and 
June 2014) to put constraints on the mass of KOI-12.01. 
The 2014 observations are of better quality than the 2011 
ones, which may be affected by a systematic shift reminis¬ 
cent of p revious S OPHIE observations in High-Efficiency 
mode (Hebrard et al.|[20l3|. Contrary to the 2014 dataset, 


Table 2. Mid-transit epoch of the LC and SC Kepler light- 
curves of KOI-12 used in our photometric analysis. 


BJD-2400000 

Short cadence Long cadence 


55033.16169 

55051.01692 

55104.58262 

55122.43786 

55140.29309 

55193.85879 

55211.71402 

55229.56926 

55283.13496 

55300.99019 

55318.84542 

55336.70066 

55390.26636 

55408.12159 

55425.97682 

55479.54252 

55497.39776 

55515.25299 

55568.81869 

55586.67392 

55604.52916 

55658.09486 

55675.95009 

55693.80532 

55747.37102 

55765.22626 

55783.08149 

55800.93672 

55836.64719 

55854.50242 

55872.35766 

55890.21289 

55943.77859 

55961.63382 

55979.48905 

56033.05475 

56050.90999 

56068.76522 

56140.18615 

56158.04139 

56211.60709 

56229.46232 

56265.17279 

56336.59372 

56354.44895 

56408.01465 


54979.59599 

54997.45122 

55068.87216 

55086.72739 

55158.14832 

55176.00356 

55247.42449 

55265.27972 

55354.55589 

55443.83206 

55461.68729 

55533.10822 

55550.96346 

55622.38439 

55711.66056 

55729.51579 

55818.79196 

55908.06812 

55925.92335 

55997.34429 

56086.62045 

56104.47569 

56175.89662 

56193.75185 

56283.02802 

56300.88325 

56372.30419 

56390.15942 


Notes: Kepler transit epochs 

are given in Barycentric Dynamical Time. 


no radial velocity standard stars were monitored during 
the 2011 observations, and we could not find any adequate 
correction for this possible shift. We thus used only obser¬ 
vations taken in 2012, 2013 and 2014. 

The transit of KOI-12.01 lasts 7.5 hours, which makes dif¬ 
ficult the observation in good conditions of the whole event 
during a night with a ground-based telescope. Run 2012 
was executed over 3 nights: the second night offers a good 
coverage of the whole transit, while reference measurements 
were taken outside the transit during the previous and the 
following nights. Run 2013 was executed over 2 nights: only 
the first part of the transit could be observed during the 
first night with partial coverage, although reference obser- 
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vations were also secured immediately before the transit, 
with an additional measurement the following night. Run 
2014 is made of two measurements taken a few days apart, 
close to the quadrature phases, to constrain the Keplerian 
semi-amplitude. Run 2013 was obtained in poor weather 
conditions and provide data of lower accuracy than Runs 
2012 and 2014. The log of the transit observations is re¬ 
ported in Table [3] 

SOPHIE is a cross-dispersed, environmentally stabilized 
echelle spectrograph de dicate d to high-precision radial ve¬ 
locity measurements (Bouchy et al. i2009). The light is dis¬ 
persed on 39 spectral orders from 387 to 694 nrn (labelled 
0 to 38 from blue to red wavelengths). SOPHIE data were 
acquired in High-Efficiency mode (resolution power A/A A 
= 40000) with exposures ranging from 600 to 1800 sec¬ 
onds depending on weather conditions, to ensure a signal- 
to-noise ratio per pixel at 550 nm as constant as possible 
(S/N 550 ~35). Spectra were extracted using the standard 
DRS SOPHIE pipeline, and passed t hrough weighted cross¬ 


correlation with a numerical mask (Baranne et al. 1996 


Pepe et al. 2002). The pipeline does not normalize the 


spectra prior to computing the cross-correlation in order to 
keep a natural weighting by the flux. The resulting cross¬ 
correlation functions (CCF) are fitted with Gaussians to 
retrieve the radial velocities (Sect. [5]), or fitted with a line- 
profile model to perform tomography (Sect. pi). We found 
that the quality of these fits was improved when removing 
some low-S/N spectral orders from the cross-correlation. As 
expected from the effects of atmospheric dispersion, blue 
spectral orders yielded the lowest S/Ns and CCFs were 
calculated using spectral orders 5-38 for Runs 2012 and 
6-38 for Run 2013. We discarded the fourth measure¬ 
ment in dataset 2013 because of its particularly low quality 
(S/N 550 =15). We tried different kinds of numerical masks 
when performing cross-correlation, and although KOI-12 
is more akin to a F-type star, the lowest dispersion on the 
order-per-order RVs was obtained using a standard G2-type 
mask. Finally, the solar light reflected by the Moon can dis¬ 
tort the shape of the CCF and shift the measured radial 
velocity. Moonlight contamination was detected in 3 expo¬ 
sures in Run 2012 and 6 exposures in Run 2013, and was 
corrected fo r by using the second fiber aperture targ eted 
on the sky (|Pollacco et al.||2008| |Hebrard et al.|[2008). We 
selected 45 SOPHIE spectra with a S/N at 550 nm greater 
than 20 to build a combined ID spectrum to be used for the 
spectroscopic analysis of the host star (Sect [3]). Individual 
exposures were set in the rest frame then co-added order 
per order using weights proportional to the S/N. The re¬ 
sulting co-added spectrum has a high S/N in the continuum 
of 150 per pixel at 550 nm (corresponding to a S/N of 548 
per resolution element). 

The SOPHIE pipeline may yield uncorrect estimates on 
radial velocities error bars when the fitted CCFs are af¬ 
fected by large rotational broadening. Since KOI-12 is a 
fast-rotator, we used instead an empirical formula well 
suited to SOPHIE radial velocity measurements, with er¬ 
ror bars defined as 3A\/FWHMccf/{CxS/'N 550 ) : where 
FWHMqcf is the full-width at half-maximum of a given 
CCF and C its contrast defined as the relative difference 
between the CCF minimum and the continuum. The radial 
velocity measurements and their uncertainties are reported 
in Table [4j Error bars on the CCFs points, needed in the 
tomography analysis, were estimated from the dispersion 


of the residuals betwe en t he CCFs and their best-fit line- 
profile model (Section [672]) . 


Table 4. Radial velocities of KOI 12. 


BJD (UTC) 

RV 

±1(7 

-2400000 

(kins -1 ) 

(kms” 1 ) 

Run 2012 56103.5518 

-21.87 

0.48 

56104.3650 

-21.26 

0.48 

56104.3832 

-20.95 

0.40 

56104.4030 

-20.88 

0.42 

56104.4201 

-21.86 

0.39 

56104.4381 

-21.96 

0.34 

56104.4580 

-22.29 

0.40 

56104.4766 

-23.00 

0.39 

56104.4884 

-22.57 

0.38 

56104.4999 

-22.76 

0.39 

56104.5110 

-21.90 

0.37 

56104.5208 

-22.11 

0.38 

56104.5299 

-22.85 

0.36 

56104.5387 

-22.98 

0.38 

56104.5480 

-22.65 

0.37 

56104.5587 

-22.80 

0.36 

56104.5675 

-22.39 

0.37 

56104.5760 

-22.76 

0.38 

56104.5861 

-22.22 

0.39 

56104.5982 

-22.38 

0.37 

56104.6155 

-22.60 

0.40 

56105.4795 

-22.26 

0.36 

56105.4886 

-22.07 

0.38 

56105.4973 

-21.99 

0.38 

56105.5351 

-22.31 

0.38 

56105.5453 

-21.99 

0.38 

56105.5568 

-22.60 

0.36 

Run 2013 56461.3690 

-24.76 

0.41 

56461.3842 

-24.57 

0.37 

56461.3992 

-24.35 

0.38 

56461.4126 

-24.72 

0.36 

56461.4259 

-25.88 

0.73 

56461.4409 

-24.35 

0.39 

56461.4559 

-23.00 

0.54 

56461.4710 

-23.16 

0.51 

56461.5313 

-24.25 

0.36 

56461.5401 

-24.12 

0.37 

56461.5488 

-24.29 

0.39 

56461.5575 

-24.02 

0.38 

56461.5662 

-24.55 

0.39 

56461.5747 

-24.24 

0.39 

56461.5831 

-23.93 

0.39 

56461.5915 

-24.07 

0.39 

56461.6010 

-24.36 

0.38 

56462.5911 

-24.24 

0.69 

Run 2014 56807.5694 

-24.46 

0.23 

56814.4686 

-24.39 

0.19 

Note: Variations in the asymmetric continuum of the 

CCFs (Fig. |4|, instrumental effects, 

or the presence 

of a massive distant companion to KOI-12 may cause 

the differences in systemic velocities between the 


three runs 


3. Spectral analysis of the host star 


Boisse et al. (2010) calibrated the relation between usim* 
and the widtn <7 of the CCF in the high-efficiency (HE) 
mode of the SOPHIE spectrograph. The value for a 
depends on the B-V magnitude (0.22±0.12 for KOI-12; 
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Table 3. Log of KOI-12.01 transit observations with SOPHIE. 


Run 

Transit mid-time (UT) 

Exposures* 

S/N 450* 

S/N.550* 

S/N 650* 

2012 

2012-06-25 at 23h28 

20 

19 

34 

35 

2013 

2013-06-18 at 01h58 

12 

18 

34 

35 


f: number of individual exposures taken during the transit. 

*: medians of the signal-to-noise ratio per pixel during the transit at 450 A, 550 A, and 650 A' 


H 0 g et al.|2000) and the full-width at half-maximum of the 
(JCEs. The average value of FWHMc^f for the two runs is 
87.5±1.9kms _1 and yields wsim*=60.0±1.5kms -1 . This 
calibration is however valid mostly for solar-type stars with 
a low or a moderate usim*. We thus checked the vsim* 


value by using the F ourier transform method (see Simon- 
Dfaz &; Herrero|2007 and reference therein). From the first 
zero in the Fourier transform of various isolated spectral 
lines in the the co-added SOPHIE spectrum, we found 
= 66.0±2.5kms ~ 1 and v macro =18. l±5.5kms _1 , 


usim* 


consistent with the value wsim* = 70.0±5.0kms 1 ob¬ 


tained by Lillo-Box et al. (2015). The final value for the 
projected rotational velocity will be given in Sect. |6.3[ 
where it is determined more accurately using Doppler 
tomography. 

We derived the star’s atmospheric parameters from 
the co-added SOPHIE spectrum o f K O I-12 usin g the 


spectral analysis package SME 2.1 (Valenti & Piskunov 


4. Photometry analysis 

The normalized Kepler light c urves were fitte d by a transit 
model using the PASTIS code ( Diaz et al.|[2014 1. The anal¬ 
ysis also includes the fit of theSED and stellar evolution 
tracks to determine cohe rent ste llar parameters. We used 
the Dartmouth (Dotter et al.j;2008) stellar evolution tracks 
as input for the stellar parameters. 

For each Kepler light curve we included the out-of-transit 
flux and the contamination factors as free parameters. We 
also account for additional sources of Gaussian noise in 
the light curves and SED by fitting a jitter value to each 


1996 Valenti & Fischer 2005). The fast rotation of the 
star prevents us to use alternative methods based on 
equivalent width measurements such as VWA. SME 
iteratively determines the fundamental stellar parameters 
by performing non-linear least-squares fit of synthetic 
spectra to the observed spectrum. The temperature was 
estimated on the extended wings of the Hydrogen Balmer 
lines. The other atmospheric parameters, logg*, the 
metallicity, but also the various velocities, usim*, v m i cro , 
and Vmacroi were estimated on a large number of metal 
lines located in different spectral windows. We found 
T eff =6800±120K, log g*=4.34±0.15, [Fe/ff]=0.09±0.15, 
usim*=66.3±1.2 kms -1 , v m i cro =2.3±0.3kms ~ 1 and 
Vmacro= 17.3± 1.4 km s _1 . We note that the values we 
obtained for the projected velocity due to rotation, 
vsi m*, and the velocity due to convection, v macro , are in 
agreement with the values we derived from the Fourier 
Transform of isolated spectral lines. From the values of 
the effective temperature, log < 7 *, and metallicity used 
as input values an d the Starevol evolutionary tracks 
Lagarde et al.|[2012 Palacios, private communication) we 
’erived the host-star’s fundamental parameters, its mass 
Af* = 1.5lioo9M Q , radius I?*=1.41±O.O6I?0, and age 
1.4±0.3 Gyr. 


dataset. This is especially appropriate for the Kepler data 
since the star is located on different CCDs each season. To 
accou nt for the wid ening of the ingress and egress dura¬ 
tions ((Kipping 2010), we modeled the LC light-curves us¬ 
ing an oversamp ling factor of 10 suitable to the Kepler data 
(Bonomo et al. 20141. We fitted the free parameters using 
a Metropolis-H asting Markov c hain Monte Carlo (MCMC) 
algor ithm (e.g. Tegmark et al.||2004 ) with an adaptive step 
size (Ford 2006| ). To better sample the posterior distribu¬ 
tion in the case of non-linear correlations between parame¬ 
ters, we applied an adaptive principal component analysis 
to the chain s and j um ped th e parameters in an uncorre¬ 
lated space (Dfaz et al.2014). For most of the parameters 
of the MCMC we used non-informative priors (uniform or 
Jeffreys distributions). Exceptions are the stellar parame¬ 
ters T e g and [Fe/H] derived from the above spectral anal¬ 
ysis, p* obtained from log p* and the evolution tracks, and 
the orbital periods and phases of the planets for which we 
used as priors the Kepler values with error bars increased 
by a factor of 100 to be conservative. We used Gaussian pri¬ 
ors for the contamination factors, centered on the MAST 


value (Sect. 2.1) and with a width of 0.2%. The priors on 


the stellar density were used to perform fits with eccentric 
orbits. 

The system was analyzed with 10 chains leading to a to¬ 
tal of 10 6 steps. Each chain was started at random points 
drawn from the joint prior. All chains converged to the same 
solution. We thinned the converged sub-chains using the 
correlation length. We finally merged the thinned chains, 
which left us with a total of more than 1000 independent 
samples of the posterior distribution. The resulting theo¬ 
retical transit light curve is displayed in Fig. [l] Because of 
an imperfect detrending or variations of the activity level 
of the star, one of the seventy-four transits was found to be 
abnormally deeper. Nonetheless, including it in the fit had 
no significant effect on the results. The best-fit parameters 
and their 68.3 % confidence intervals are listed in Table [5J 
We did not detect a significant eccentricity, with a 3er upper 
limit of 0.72. Other parameter values were found to be in 
agreement with those derived for a circular fit, performed 
on the transit parameters without adjusting the stellar pa¬ 
rameters (Table [ 5 ]) . Hereafter we choose parameters from 
the circular fits as final values. The results obtained in this 
section are also taken as final values for the stellar param¬ 
eters. 

5. Radial velocity analysis 

5.1. Keplerian fit 

We fitted radial velocities outside of the transit in Runs 
2012, 2013 and 2014 simultaneously, using a Keplerian or¬ 
bit. The orbital period P and transit epoch To were fixed 
to their photometric values (Table [H]). We computed the % 2 
of the fit on a grid scanning a large range of values for the 
semi-amplitude K of the radial velocity variations, and the 
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Table 5. Transit light curve and radial velocity analysis. 


Parameter 

Symbol 

Value 

Unit 

Eccentric orbit 

Orbital parameters 


17.8552333i0.9x 10” 6 


Orbital period 

P 

day 

Transit epoch 

T 0 

2454979.59599±2.7x 10” 4 

BJD 

Eccentricity 

e 

<0.72 (3a) 


Argument of periastron 

UJ 

228±120 

deg 

Scaled semi-major axis 

dp / R* 

20 .Oil.5 


Semi-major axis 

dp 

0.151i3x 10” 3 

au 

Orbital inclination 

ip 

88.95±o:25 

deg 

Impact parameter 

b 

0.3629±6.4x 10” 3 


Planet-to-star radii ratio 

Rp/R * 

0.09049i8xl0” 5 


Planet radius 

R p 

1.43F0.13 

R'jup 

Transit duration 

td 

a 7 oo + 1-8x10- 2 

D. 1 OO-0.7X10- 2 

hours 

Stellar reflex velocity 

K 

<0.8(3cr) 

km s” 1 

Mass of KOI-12.01 

M p 

< 10 (3<t) 

Mj up 

Stellar parameters 




Density 

P * 

0.335 i 0.080 

Psun 

Radius 

R* 

1.63 i 0.15 

Rsun 

Age 

taUi, 

1.5 i 0.5 

Gyr 

Mass 

M* 

1.452 i 0.093 

M-sun 

Distance 

D 

426 i 40 

pc 

Effective temperature 

T e ff 

6820 ± 120 

K 

Linrb-darkening coefficients 

€ a 

0.334F0.007 



Cb 

0.124F0.013 


Circular orbit 

Orbital parameters 




Orbital period 

P 

17.8552332il.0xl0” 6 

day 

Transit epoch 

To 

2454979.59601i0.5x 10” 4 

BJD 

Scaled semi-major axis 

dp/ R-k 

18.84i0.04 


Semi-major axis 

dp 

0.151i3x 10” 3 

au 

Orbital inclination 

ip 

88.90i0.02 

deg 

Impact parameter 

b 

0.362i7x 10” 3 


Planet-to-star radii ratio 

Rp / R* 

0.09049i8x 10” 5 


Planet radius 

R P 

1.44i0.13 

R'./itp 

Transit duration 

td 

6.7827i2.6x 10” 3 

hours 

Stellar reflex velocity 

K 

< 0.51 (3ct) 

km s” 1 

Mass of KOI-12.01 

M p 

< 8.7 (3er) 

Mjup 

Stellar parameters 




Limb-darkening coefficients 

€ a 

0.335i0.008 



Sb 

0.123i0.014 



Notes: All parameters are derived from the analysis of the transit light curve (Sect. |4h , except for the limits 


on the Keplerian semi-amplitude and the mass of KOI-12.01 derived from the radialvelocity analysis (Sect. 5.1l. 
Notes: The derived Kepler transit epochs are given in Barycentric Dynamical Time. 


systemic radial velocities for each dataset 72012 , 72013 , 72014 - 
We adjusted independently the systemic velocities because 
of possible variations in the asymmetric continuum of the 
CCFs, or resid ual instrumental drifts, between the three 
runs (Sect. 2.2). Once the minimum \ 2 and corresponding 


best values for these parameters were obtained, we calcu¬ 
lated their error bars from an analysis of the y 2 variation 


(see e.g. Hebrard et al. 2002). 

Radial velocity measurements and their Keplerian model 
for a circular orbit are plotted in Fig. [2] We found sys¬ 
temic radial velocities of -22.17±0.09, -24.57±0.10 and - 
24.44i0.09kms” 1 for Runs 2012, 2013 and 2014 respec¬ 
tively. By comparison th e RV measured in 20 11 are in 
the order of -18.5 kms” 1 (Santerne et al. 2012a). The sig¬ 
nificant variations in systemic velocity between the runs 
may be caused by variations in the asymmetric continuum 


of the CCFs, residual instrumental drifts, or the presence 
of a distant massive companion to KOI-12. However the 
small number of measurements over long timescales, the 
lack of RV standard stars in 201 1 and possible systematics 
in SOPHIE HE mode (Sect. 2.2) prevent us from conclud¬ 
ing. The effects of the asymmetric continuum were ho weve r 
taken into account in the tomography analysis (Sect |6~3] ), 
allowing us to refine the value for the systemic velocity of 
KOI-12 in Run 2012. The precision of our RV measure¬ 
ments and the sampling of the orbit are insufficient for 
a significant detection of the star’s reflex motion, and we 
derived a 3cr upper limit of 0.51 kms” 1 on the Keplerian 
semi-amplitude K (Table [5]). This parameter is mainly con¬ 
strained by measurements in Run 2014, which are closer 
to quadrature-phase than measurements in Runs 2012 and 
2013 taken near the transit epoch. Using the orbital in- 
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clination z p and stellar mass M* obtained in Sect. [4j to¬ 
gether with their uncertainties, the limit for K corresponds 
to a 3 a upper mass limit M p =8.7 Mj up for KOI-12.01. This 
val ue is more stringent tha n the limit of 26.7 Mj up obtained 
by Santerne et ah (2012aI, as the more recent datasets we 
used in Sect. 2.2 are oFTTetter quality with more measure¬ 
ments. It is also consistent with the maximum proj ected 
mass of 25.2 Mj up determined by Lillo-Box et al. (20151. 


We performed the same analysis assuming an eccentric or¬ 
bit, setting the eccentricity to its 3er upper limit of 0.72. 
We found that the 3 ct upper limits on K and M p increased 
to 0.8 kms -1 and 10Mj up , which is taken hereafter as the 
final upper limit on the mass of KOI-12.01. 


5.2. Rossiter-McLaughlin anomaly 



- 0.06 - 0.04 - 0.02 0.00 0.02 0.04 0.06 

Transit phase 


Spectra in Run 2013 were acquired in poor w eather 
conditions, and their tomographic analysis (see Sect. 6.4.2) 
allowed systematic effects to be identified in two of the 
radial velocity measurements (see Fig. [2]). In addition, most 
of Run 2013 measurements are grouped near the transit 
center, which makes it difficult to constrain the shape 
of the RM anomaly. In contrast, Run 2012 has a good 
sampling of the whole transit with data of good quality. We 
thus applied the analytical approach developed by Ohta| 


et al. (2005) to model the form of the RM anomaly in this 
run only. The semi-amplitude was fixed to Okms' 1 and 
the syst emic velocity was fixed to its best-fit value derived 
in Sect. [oTT| for a circular orbit. The stellar limb-darkening 
coefficients, the transit parameters R v /R+, a p /i?* and ip, 
P and To were fixed to their photometric values (Table pi). 
The best-fit values and corresponding error bars for the 
sky-projected obliquity A and the projected stellar rotation 
velocity risinz* were determined in the same way as in 


Sect. |5.1| using analysis. 

The best model for the RM anomaly is shown in Fig. [2] 
It was obtained with usinz*=117±25kms _1 , which is 2er 
higher than the value derived from the spectral analysis in 
Sect. |3l The analytical formulae developed by |Ohta et al.| 
( 2005p do not provide a good estimate of the projected 
stellar r otation velocity f or rapidly rotating stars like 
KOI-12 ( Hirano et al.|20fo|). Yet, we also analyzed the RM 
anomaly using the Gimenez (2006) equations and found a 
similar value with z;smz*=107±23 kms -1 . Both methods do 
not take into account the effects of macroturbulence, which 
are expected to be important for KOI-12 (Sect|3]). This 
may have biased the estimates of the projected rotational 
velocity toward higher values, although a macroturbulence 
in the order of 18kms -1 cannot explain an overestimation 
by about 50 kms -1 . The value for usinz* will be determined 
more accurately using the technique of Doppler tomogra¬ 
phy in Sect. [6] Nonetheless, anal ysis of the RM ano maly 


using the analytical formulae of Ohta et al. (2005) and 


Gimenez (20061) provides a useful first-order estimate of the 


obliquity value, and we measured A = 5.8lti5.'5 using both 
methods. The fitted obliquity remained consistent within 
its uncertainties when we increased to their 3er upper 


limits the Keplerian semi-amplitude K (A = -14.7li5;s°) 

or the eccentricity e (A = lTlil'i °), as in both cases the 
Keplerian velocity variations remain small on the transit 
duration timescale (see Fig. [2]). Finally, we checked the 
consistency of Run 2013 with the RM anomaly detected 
in Run 2012 by fitting them simultaneously, assuming no 


Fig. 3. Upper panel: Error-weighted average of Runs 2012 
and 2013 RV measurements (black points), following the 
subtraction of the systemic velocities and assuming no re¬ 
flex stellar motion (K =0 km s -1 ). The RM adj ustment with 

A=16. 0 ^ 10.4 and t>sinz*=123.5±22.5kms -1 is overplotted 
with a red line. Keplerian residuals are first calculated in 
each run and then grouped in common phase bins. We 
excluded the two anomalous RV measurements in dataset 
2013 (Fig. [2]). Lower panel: Residuals from the Keplerian + 
RM fit, without binning. Green points correspond to Run 
2012, orange points to Run 2013. 


stellar reflex motion. The best RM model is plotted in 

a _i_9 g 0 

and corresponds to an obliquity of 16.0_ 10.4 and 
a projected rotational velocity of 123.5±22.5kms -1 , in 
agreement with the values obtained using Run 2012 only. 
The two RV measurements affected by systematic effects 
in dataset 2013 were excluded from the simultaneous fit, as 
they are strongly redshifted and force the stellar rotational 
velocity to a higher value (see the tomographic analysis 
and Fig. [8]). Nonetheless, we found that including them in 
the fit does not change significantly the derived obliquity. 
We thus conclude to the detection of the RM anomaly 
of KOI-12.01, which shows a prograde orbit consistent 
with an aligned system in both transit observations. We 
note, though, that the large uncertainties derived for the 
obliquity do not rule out a small misalignment. 


6. Doppler Tomography 

6.1. Parameter fitting method 

We then used Doppler tomography to obtain complemen¬ 
tary values for the system properties and investigate their 
consistency with the fit to the RM anomaly. Although this 
latter analysis showed a prograde orbit and excluded a high 
misalignment, it could not indeed place strong constraints 
on the obliquity and yielded an anomalously high value 
for the projected stellar rotation velocity. We expect a 
more reliable value for this parameter to be derived using 
tomography, as this method is well suited to fast rotators 


such as KOI-12 ( 

Collier Cameron et al. 

2010b 

) and takes 

into account the effects of macroturbu 

ence through the 


direct modeling of the local line profile. We used prior 
constraints from photometry to increase the accuracy of 
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Transit phase 


Fig. 2 . Radial velocity measurements during Runs 2012, 2013 and 2014 (black points, from top to bottom), and their 
circular Keplerian + RM fit for a null semi-amplitude (red line) and its 3cr upper limit (blue line). The Keplerian fit was 
performed on all measurements outside of the transit, and yields a total dispersion on the residuals of 0.21 and 0.28 km s -1 , 
respectively. Variations in the asymmetric continuum of the CCFs (Fig. [4]) may cause the apparent differences in systemic 
velocities between the three runs. The RM model was adjusted to dataset 2012 only, but is consistent with dataset 2013. 
The two anomalous measurements in Run 2013 are the most redshifted RV values near ingress (vertical dashed lines 
show the times of mid-transit, first, and fourth contacts). 


our results. For the reasons given in Sect. |5.2| we based 
our analysis on the spectra secured during Run 2012, and 
in a second time we checke d the consistency of dataset 
2013 with our results (Sect. 6.4.2). The CCFs obtained 
with the SOPHIE pipelin e were fitted u sing the technique 
developped by Collier Cameron et al. (2010a), which 
relies on the decomposition of the CCF into its different 
components: the stellar line profile model outside the 
transit is a limb-darkened rotation profile convolved with 
a Gaussian corresponding to the intrinsic photospheric 
line profile plus instrumental broadening. The planet 
occultation is modeled as a Gaussian ’bump’ in the stellar 
line profile, whose spectral location depends on the planet 
position in front of the stellar disk during the transit. The 
continuum of KOI-12 CCFs is asymmetric and strongly 
tilted (Fig. [4]). While stellar activity, moonlight and 
telluric contamination, and random alignments between 
photospheric lines and mask lines at arbitrary RV shifts 
are known to create anomalies in the profile of the CCF 
(the so-called sidelobe patterns), here the observed tilt 
may be caused by an imperfect removal of the instrument 
blaze function. While this pattern may decrease the quality 
of the CCFs adjustment, we noted that it remains fixed 


over a few consecutive nights and can thus be corrected for 
by subtracting the average of the diff erence between each 
CC F and its model line profile (Collier Cameron et al. 
2010a). To prevent removing inadvertently KO1-12.01 


signature, however, only observations outside the transit 
were used to calculate this correction. 

The tomographic model depends on the same pa¬ 
rameters as the Keplerian + RM model (Sect. 5l), 
with the addition of the local line profile width 
s=FWHM/(2 \j2ln 2 vsini*), expr essed in uni ts of the p ro- 
jected stellar rotational velocity (Collier Cameron et al. 
2010a). FWHM is the full-width aFhalf-maximum of the 


non-rotating stellar line profile convolved with the instru¬ 
ment profilcQ The effects of macroturbulent velocity on the 
width of the local line are thus already accounted for by 
the FWHM, independently of usirn*. We fixed the period 
and the mid-transit time to their values determined with a 
high precision using the Kepler photomet ry ( Sect.j4|). The 
semi-amplitude was set to Okms -1 (Sect. 5.1) and we used 
a linear limb-darkening law with coefficient e. We discuss 


The FWHMccr introduced in Sect. 2.2 is instead the full- 


width at half-maximum of the rotationally-broadened CCF 
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Fig. 4. Mean of the CCFs observed outside of the transit in 
Run 2012, as a function of radial velocities relative to the 
star. Before calculating the average, the CCFs in each run 
were shifted from the velocity frame of the Solar system 
barycenter where they are computed, to the frame of KOI- 
12 using the best-fit systemic velocities derived in Sect. |5.1| 
Note the strong tilt of the asymmetric continuum. 


the influence of these parameters values and assumptions 
on the results in Sect. All other parameters were ad¬ 

justed using the following merit function: 


X 


nccF riv r 

= ££ 


fi,j(obs ) - (mod) 


1 2 


+ 


£ 

& p / R* j^p jRp / R* 


%tomo •Kphoto 


Z'photo 


( 1 ) 


where fij is the flux at velocity point j in the ith observed 
or model CCFs. The error on the flux, <7j, is supposed con¬ 
stant for a given CCF. Photometry cumulated over seventy- 
four Kepler transits is expected to provide a more precise 
determination for the radius, semi-major axis and orbital 
inclination of KOI-12.01 than tomography. The model value 
of these parameters Xt omo is constrained directly through 
the y 1 2 3 using the best-fit value Xph oto ± cr Xphoto from the 
transit light-curve analysis (Table pi) . The fit to the data is 
performed in three steps: 


1. To identify all possible local minima in the parameters 
space, we analyze y 2 variations. A given parameter is 
pegged at various trial values, and for each trial value 
we run a Levenberg-Marquardt minimization algorithm, 
allowing all other parameters to vary freely. 

2. We use the residuals between the CCFs and their best- 
fit model profile from step 1 to estimate the noise on the 


CCFs. The procedure is described in details in Sect. |6.2 
Note that the value of cr* has little influence on the re¬ 
sults of step 1. 

3. We refine the best-fit parameters values and calculate 
their uncertainties using MCMC. The Markov chain is 
constructed by repeatedly computing the merit function 
for a sequence of parameter values. At each step of the 
chain, the difference Ay 2 is calculated and the new set 
of parameters is accepted if Ay 2 <0, or if Ay 2 >0 with a 
probability equal to the ratio of likelihoods. The next set 
of parameters is computed by applying Gaussian pertur¬ 
bations to the preceding set: x i+ i = aq + Gauss(0, a x ). 
The standard deviation a x is specific to each parame¬ 
ter, and is reevaluated for all parameters with a con¬ 
stant frequency, so that a change Xj+i = aq + a x would 
correspond to Ay 2 = 1. This ensures that the space 
of each parameter is correctly sampled and that their 


posterior probability distributions have the correct vari¬ 
ances. We start several independent chains from differ¬ 
ent sets of parameters chosen randomly near the local 
minima found in step 1. This limits the burn-in phase 
to only a few hundred points, after which the chain is in 
a steady state and all accepted steps are used to com¬ 
pute the probability distributions. Although the chains 
are let free to explore all values for the orbital incli¬ 
nation and obliquity, we use the following conventions 
to compute the final distributions: 0° < i p < 90° ; - 
180° < A <180°. Accordingly, all couples ( i p >90°, A) 
are tranformed into (180° — i p , —A) as the two possi¬ 
bilities are equivalent with respect to the tomography 
model equations. The medians of the distributions are 
chosen as the final best-fit values for the model param¬ 
eters. Their ler uncertainties are obtained by finding 
the intervals on both sides of the median that contain 
34.15% of the accepted steps. 


6.2. Determination of CCFs errors 


Because errors on the CCF data points are unknown, we 
attributed at first to each pixel of a given CCF a con¬ 
stant error equal to the dispersion of the residuals between 
this CCF and its best-fit model profile. However the CCFs 
were calculated with the SOPHIE pipeline at a velocity 
resolution of 0.5kms _1 , while the spectra were observed 
in the HE mode at an instrumental resolution of about 
7.5kms _1 . The residuals were thus found to be strongly 
correlated, which can lead to underestimation of error bars 
on the derived parameters. To derive reliable uncertainties 
on the model parameters using y 2 statistics, errors on the 
CCF pixels must be independent and follow a Gaussian 
random distribution (a ’white’ noise). We thus decided to 
retrieve directly the non-correlated gaussian component of 
the CCFs noi se usi n g an analysis of the residuals variance 
inspired from Pont et al. (2006): 

1. We calculate the CCF residuals defined as the difference 
between a CCF and its best-fit mod el obtained in step 
1 of the fitting procedure (Sect. 6.1). 

2. We average the residuals within an interval of length 
ribin pixels (the binning factor). This interval slides 
along the entire residual range, and we record the resid¬ 
ual means for every possible position where the interval 
contains a different set of 'ribin successive pixels. 

3. To each value of the binning factor corresponds a dis¬ 
tribution of residuals means, for which we calculate the 
standard deviation a{nbi n )- This parameter represents 
the characteristic dispersion of the CCF pixels binned 
at a given resolution. This dispersion is the result of 
the combination of the intrinsic Gaussian noise in indi¬ 
vidual measurements and the correlation between them. 


After calculating the standard deviation for every CCF, we 
found empirically that it is always well represented by a 
quartic harmonic combination of a white and a red noise 
components: 


rr {ribin) 


Rbin | 
®Uncorr J 


2 


+ 


1 

2 

®Corr 


1 

2 


( 2 ) 


where <Jjjncorr/ V ribin can be understood as the result of 
the intrinsic uncorrelated noise after the binning of ribin 
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Fig. 5. Standard deviation of velocity-binned residuals be¬ 
tween a CCF and its best-fit model, as a function of the bin¬ 
ning factor (black diamonds). A single pixel is 0.5kms _1 
wide. The best fit (black line) is a harmonic combination 
of an uncorrelated noise component, inversely proportional 
to the square root of the binning factor (blue line), and a 
constant correlated noise component (red line). 


pixels, and crcorr is a constant term characterizing the 
correlation between the binned pixels. A typical curve for 
the standard deviation as a function of the binning factor is 
shown in Fig[bJ To validate Eq. [2] we created thousands of 
synthetic CCF residual tables in which each point is drawn 
from a Gaussian random distribution. These residuals were 
interpolated onto a velocity table with higher resolution 
to make them equivalent to the CCF residuals resulting 
from the model fits to the observed CCFs, by creating 
correlation between adjacent pixels as observed in the 
data. For each synthetic table we then calculated cr(jibin) 
using the above method, and verified that it was always 
best described by Eq. [2j 

As can be seen in Fig [5] the uncorrelated noise in the 
CCFs becomes dominant for binning factors higher than 
about 15, because the velocity width of the bins is then 
larger than SOPHIE instrumental resolution. In previous 
tomography studies, the issue of correlated noise was 
circumvente d by binning CCF s at the instrumental reso¬ 
lution (e.g. Miller et al. 2010). To avoid possible loss in 
resolution, we did not bin the data and attributed to all 
pixels in a given CCF the same uncorrelated noise crjj ncorr . 
This value represents the theoretical intrinsic Gaussian 
noise (the blue line in Fig pH) of an unbinned pixel. With 
CCFs errors now corresponding to uncorrelated noise only, 
uncertainties on the model parameters derived from the fit 
of the CCFs can be correctly estimated using y 2 statistics, 
independently of the CCF sampling. 


6.3. Results 

In this section we present our results from the tomographic 
analysis of Run 2012 time-series CCFs, which are plotted in 
the top panel of Fig. [6j Much information is readily avail¬ 
able in the residual maps obtained by subtraction of the 
out-of transit model stellar line profile (Fig. fbj upper mid¬ 
dle panel). The transit is revealed as a bright streak con¬ 
fined between ± usini*, and whose spectral positions over 
time correspond to the radial velocity of the stellar sur¬ 


face regions occulted by KOI-12.01. For comparison, Fig. [6] 
(lower middle panel) shows the signature of a transiting 
planet that appears as a straight line in the residual maps. 
Indeed, its velocity can be written as v p = vsini * x±, where 
x± is the distance between the region occulted at orbital 
phase (j) and the stellar spin axis, 

x± = ap/Ri, (cos(A) sin(2 7T(^) — sin(A) cos(i p ) cos(2 7r 0)). 

( 3 ) 

During the transit, </> ~ 0 and the velocity of the signature 
approximates to a linear function of the orbital phase 


v p = vsini * dp/Ri, (cos(A) 2 7r <f> — sin(A) cos(i p )). (4) 


The signature of KOI-12.01 travels from negative to 
positive velocities during the transit, which shows unam¬ 
biguously that it is on a prograde orbit, occulting first the 
blueshifted and then the redshifted regions of the stellar 
disk. As can be seen in Fig. [6] (upper middle panel), the 
signature trajectory is roughly symmetric, ruling out a high 
misalignement consistently with the analysis of the RM 
anomaly in Sect. |5.2| However there is a hint that ingress 
occurs at a slightly higher absolute velocity (^60kms _1 ) 
than egress (~50kms _1 ), the former being close to the 
best-fit projected stellar rotation velocity (Table [6]). This 
would be consistent with a moderately inclined orbit, the 
transit beginning near the equatorial plane and ending 
at a higher latitude. Finally, faint features are visible in 
the overall residual maps after subtraction of the model 
planet signature (bo ttom pan el in Fig, [dl). As in the 
case of WASP-33 (Collier Cameron "etaT 2010b), such 
features could be due to oscillations of the stellar surface 
which move alternatively inward or outward depending 
on the longitude. However, the prograde motion of these 
signatures and their linearity in the phase/velocity residual 
maps suggest the presence of active regions carried around 
by the rotation of the star. 


Only one set of parameters values was found to 
reproduce well the ob serv ations during step 1 of the 
fitting procedure (Sect. 6.1). We ran eight MCMC chains 
of 10 4 accepted steps, with an average acceptance rate 
of 25%. The posterior probability distributions for all 
model parameters are shown in Fig. [TJ along with the 
marginalized ID distributions which are well represented 
by Gaussians and allowed us to derive tight constraints 
on the inferred best-fit values (Table. [6]). The combined 
analysis and the use of priors on a p /R*, i? p /i?*, and 
i p reduced the uncertainties on these parameters. We 

found usim*=60.0^o.'g kms~ 4 , lower by about 2 a than 
the value obtained using the Fourier transform method 
but in good agreement with the value derived from the 


Boisse et al. 2010 relationship (Sect, 
that KOl-12 


This confirms 
rapidly rotating star 


is one of the most 
known to host a transiting planet so far. With such a 
large rotational broadening, the value derived for the local 
line width s=0.076±5x 10 -3 corresponds to a FWHM for 
the local stellar and instrumental profile of about 18% of 
the proje cted rotation velocity (see the definition of s in 
Sect. 16.1 ). 


i 2 o 

The best-fit value for the obliquity A=12.5_2.g is consistent 
with the RM anomaly fit, and the visual analysis of the 
residual maps that hinted to an asymmetry in the signature 
of KOI-12.01 (see Fig. [6]). The reduced uncertainties allow 
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Velocity (km/s) 


Fig. 6. Maps of the time-series CCFs from Run 2012, as a function of radial velocity relative to the star (in abscissa) and 
orbital phase (in ordinate, increasing vertically). Flux values increase from dark red to white. Vertical dashed white lines 
are plotted at ± vsi m*, and white diamonds indicate the time of the 1 st and 4 tft contacts. Top panel : CCFs produc ed b y 
the SOPHIE pipeline. Dashed white areas separate the different datasets taken over three successive nights (see Sect. 2.2). 
The tilt of the continuum is visible as a flux increase from negative to positive velocities (see also Fig. [4]), and does not 
vary significantly from one night to the other. Upper middle panel: Map of the transit residuals after subtracting the 
model stellar profile. The signature of KOI-12.01 is the bright and wide feature that crosses the entire width of the 
line profile from ingress to egress. Lower middle panel: Best-fit model for KOI-12.01 transiting signature, emphasized 
by the dark brown line (extended beyond the first/fourth contact and ± i>sini* for the sake of clarity) obtained with A 

= 12.5_2’9- This is consistent with the asymmetry of the signature visible in the above panel, whereas an aligned orbit 
would correspond to the lighter brown line. Bottom panel: Overall residual map after the further subtraction of the model 
planet signature. The faint, alternating bright and dark parallel streaks may be caused by stellar intrinsic variability. 


us to conclude to the low-obliquity of KOI-12.01 orbit. 

The systemic veloci ty o btained from Gaussian-fitted RV 
measurements (Sect. 5.1) is significantly different from the 
value derived from tomography. This is likely because the 
latter method is based on a direct modeling of the CCFs 
which takes into account the asymmetric non-flat contin¬ 
uum of the CCF (Fig. [4]), while a Gaussian adjustment 
with a flat continuum model introduces systematic shifts 
in the estimation of the velocity. Nonetheless this does not 
prevent to correctly sample the Keplerian velocity curve, 
as the patterns remain similar over a few nights and the 
velocity shifts are thus constant in a given dataset. 


Finally, the diagrams in Fig. [7] show no parameter cor¬ 
relations except between A, 72012 and vsim*, and we there¬ 


fore studied the influence of these three parameters on the 
model. The quality of the observed CCFs and the lower 
phase coverage during the first half of the transit allow for 
some indetermination in the spectral position of the plane¬ 
tary signature in the stellar reference frame. This velocity 
reference depends directly on the systemic velocity, and a 
higher value for 72012 is equivalent in the residual maps 
to a shift of the planetary trace toward negative velocities 
relative to the star. This requires in turn that the planet oc¬ 
cults the blue regions of the stellar disk for a larger amount 
of time, which can only be obtained with a higher obliq¬ 
uity since the impact parameter is tightly constrained by 
the photometry priors. There is thus a strong correlation 
between the systemic velocity and the obliquity, and an 
increase in these parameters also requires an increase in 
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Table 6. Doppler tomography analysis, with and without constraints from the photometry analysis. Pertinent values 
from Table [5] have been reported in the first column for a more direct comparison. The final values for the system 
properties are derived from the tomography + photometry analysis. 


Parameter 

Symbol 

Photometry 

Tomography 

Tomography + Photometry 

Unit 

Scaled semi-major axis 

&p / R-k 

18.84±0.04 

17.7-0.7 

18.83±0.03 f 


Orbital inclination 

ip 

88.90i0.02 

89.3io.-7 

88.902±0.015 f 

deg 

Impact parameter 

b 

0.362±7xl0“ 3 

0.2±g;? 

0.361±5xl0 -3 


Planet-to-star radii ratio 

Rp/Ri, 

0.09049±8x 10 -5 

0.099±4x 10 -3 

0.09050±6x 10~ 5 f 


Stellar rotation velocity 

vsi m* 

- 

62.9T1.4 

60.0t o °k 9 

km s" 

Local line width 

s ft 

- 

0.084±6xl0 -3 

0.076±5xl0~ 3 


Systemic velocity 

72012 

- 

-16.1±1.1 

-14.7±1.1 

km s" 1 

Sky-projected obliquity 

A 

- 

12.0t6 3 6 4 

1 0 c - ! -3,0 

12.0 -2.9 

deg 


f: These parameters are constrained with priors from photometry. 


tt : s=FWHM/(2 y/2ln2 vsini*) 


the rotational velocity since the shifted planetary signature 
must remain within the range of ± r>sini* during the transit. 


6.4. Consistency checks 

6.4.1. Model parameters 


The shape of the stellar line profile and the amplitude of the 
planetary signature near ingress and egress depend on the 
effects of limb-darkening. In the present case however, the 
quality of the data and the phase coverage/resolution do 
not allow for the determination of the limb-darkening co¬ 
efficient. All model parameters values derived in Sect. |6.3| 
thus remained within their lcr uncertainties when e was var¬ 
ied between 0.1 and 0.5. Because of the quality of the data, 
using a Voigt profile instead of a Gaussian profile to model 
the local line profile had no influence on the fit. 

We ran a tomographic fit for an eccentric orbit with e set to 
its 3 a upper limit of 0.72 and found that all model parame¬ 
ters, in particular the obliquity, remained consistent within 
their lcr uncertainties with their best-fit values derived in 
Sect. 16.31 Similar results were obtained for a circular or¬ 
bit with the Keplerian se mi-a mplitude set to its 3cr upper 
limit of 0.51 kms -1 (Sect. 5.1). Both eccentricity and semi- 
amplitude affect the shape of the Keplerian velocity curve, 
shifting the location of the stellar rest velocity in each of the 
observed CCF. This can in theory modify the best-fit val¬ 
ues for the obliquity, systemic velocit y, an d stellar rotation 
velocity for the reasons given in Sect. |6.3| Here though, RV 
variations at the time of the observations remained small 
enough to have no significant effect on the results. 


6.4.2. Run 2013 


We performed the tomographic analysis of dataset 2013 
alone to check its consistency with dataset 2012, but found 
we could not obtain a good fit to the data. Instead, we cal¬ 
culated the residual maps for R un 2 013 using the best-fit 
stellar line profile derived in Sect. 6.3 for Run 2012 (Fig.|8|. 
The transiting signature of KOI-12.01 was revealed in the 
residuals (upper panel), and is consistent with the best-fit 
planet signature from Run 2012 analysis (middle panel). 
The fit to dataset 2013 may have failed for several reasons. 
Tomography is most efficient when the entire planet tra¬ 
jectory can be sampled in the time-series CCFs, as was the 


case in Run 2012. Unfortunately, Run 2013 has a poor phase 
coverage with most measurements taken near the center of 
the transit and nearly no observations during its second 
half (see Fig. |2|). In addition, the lack of reference spec¬ 
tra outside of the transit likely hindered the correction of 
the asymmetric continuum. The quality of the adjustment 
may also have been decreased by anomalies such as the 
low fluxes in the red wing of two CCFs near ingress (see 
Fig. [8]) . These features are respo nsibl e for the spurious RV 
measurements mentioned in Sect. |5.2| since an artificial dip 
in the red wing of a CCF forces its Gaussian adjustement 
toward a higher fitted RV. Residual maps obtained with 
tomography are a powerful way to identify such anomalies, 
especially when their influence on the RV measurements 
cannot be identified directly, as was the case in Fig. [2] 


6.4.3. Tomography-only analysis 


We investigated whether tomography alone could yield 
results consistent with the photometry and photome- 
try+tomography analyses. We used the same procedure but 
removed the constraints from the pho tome tric priors. The 
first step of the fitting procedure (Sect. 6.1) showed some of 


the model parameters to yield a good fit to the data over a 
wide range of values, and we thus started the MCMC chains 
from diverse sets of acceptable parameter values. All chains 
nonetheless converged toward similar distributions, whose 
cumulated posterior probability distributions are displayed 
in Fig. [9] The resulting best-fit parameter values are given 
in Table. [6] 

We sought to understand why some of the parameters 
present strong correlations and much larger uncertainties 
than in the photometry+tomography analysis. A band of 
stellar surface parallel to the spin axis is characterized by 
a unique radial velocity v p proportional to its distance x± 
to the stellar rotation axis (Eq.E|. Neglecting the effects of 
limb-darkening, the signature ofKOI-12 at a given velocity 
v p in the CCFs can thus come from the occultation of any 
part of the stellar band at x±. As many combinations of 
a p /I?*, i p and A yield the same value of x±_ (Eq. |5]0 the 
absence of photometry priors for these parameters result in 
larger uncertainties. The degeneracy is limited by the good 


5 a p /R t and i p can be constrained independently thanks to 
the phase-dependent term in Eq. [3] 
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ap/R. i p Rp/R. vsini. s X y A 

Fig. 7. Correlation diagrams for the probability distributions of the tomographic model parameters. Priors from photom¬ 
etry were used to constrain a p /f?*, i p and The green and blue lines show the 1 and 2 a simultaneous confidence 

regions that contain respectively 39.3% and 86.5% of the accepted steps. ID histograms correspond to the distributions 
projected on the space of each line parameter. The red line and white point show median values. Units are the same as 
in Table [Gj 


phase coverage of the 2012 transit and corresponding good 
sampling of the x± values. Although the distributions for 
a p /i?*, i p and A may seem bimodal (Fig. [9]) we found that 
the two apparent modes are part of a unique region of min¬ 
imal x 2 - There is thus only one set of acceptable values for 
all parameters, albeit with large uncertainties. 

We then compared the consistency of these results with 
those obtained in the previous analyses. Best-fit values for 
Rp/Ri, and wsini* are higher by about 2er with respect 
to their values derived from the photometry+tomography 
analysis. It is unclear why we obtain this difference for 
the projected stellar rotation velocity, since it should be 
well constrained by the shape of the unocculted stellar line 
profile and by the spectral trajectory of KOI-12.01 in the 
residual map. The large value for I? p /I?*, however, may 


have several origins. First, macro-turbulence elongates the 


wings of the local stellar line profile (e.g., Hirano et al. 


2011), affecting the width of the planet bump. Because 


we assumed this profile to be gaussian, and KOI-12 has a 
large macro-turbulence velocity ( v macro ~18kms , higher 
than the resolving power of the SOPHIE spectrograph at 
7.5 km s^ 1 ), the fit may have been biased toward a deeper 
gaussian profile, ie higher values for i? p /l?*. Secondly, the 
width of the planet bump depends on photospheric turbu¬ 
lence, and can be affected by blur of the planet motion if 
the exposure time is too long. We found that four in-transit 
exposures, in the first half of Run 2012 transit, have expo¬ 
sure durations larger than about 20mn, which is the limit 
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Fig. 8. Upper panel: Residual map of dataset 2013 after 
subtraction of Run 2012 best model stellar profile, plotted 
on the same scales as in Fig. [6] for a better comparison. 
CCFs are plotted as a function of radial velocity relative 
to the star, with vertical dashed white lines at ± esini*. 
Orbital phase increases vertically, with white diamonds 
plotted at ingress and egress. The signature of KOI-12.01 is 
mainly visible as a moderately bright streak near the cen¬ 
ter of the transit, and is consistent in terms of phase and 
velocity with the best-fit model signature from Run 2012, 
displayed in the middle panel at the phases of Run 2013 
spectra. Lower panel: Residual map after the further sub¬ 
traction of Run 2012 model planet signature. Data is noisier 
at the beginning of the run, with two CCFs in particular 
(phases -0.0071 and -0.0062) showing sharp anomalous dips 
near 25kms _1 (highlighted with a green ellipse). 


we estimated for motion blur to become significant^] It is 
unclear, however, if such a small number of blurred expo¬ 
sures would be enough to bias the estimation of R p /R*. A 
last possibility may be the merging of KOI-12.01 signature 
with other features near phase 0.004 (see middle panel in 
Fig. [6]). The resulting brighter region is located in the best- 
sampled part of the time-series spectra which constrains 
much of the fit, and may thus force the amplitude of the 
model planet signature, proportional to the planet-to-star 
surface ratio, toward larger value. As for the analysis of 
Run 2013, this illustrates how tomography depends on the 
phase sampling and the quality of the spectra (e.g. statisti¬ 
cal noise, or signatures due to stellar activity), two points 
which can limit the ability of this technique to correctly 
retrieve the planet signature and the corresponding system 
properties. Except for R v /R+ and usim* the tomographic 
analysis alone, performed on a single transit, nonetheless 
yields model parameter values consistent within their 1 <j 


6 A rough estimate of the required exposure duration can be 
obtained by assuming the bump should not move by half of the 
local line profile width during the exposure 


uncertainties with the results from photomometry, its com¬ 
bination with tomography, and with the analysis of the RV 
anomaly for the obliquity (Table. [6]. 

7. The planetary nature of KOI-12.01 

To confirm the planetary nature of KOI-12.01, we consider 
here the different false positive scenarios that could have 
mimicked the photometric light curve observed with Kepler: 

1. KOI-12 and KOI-12.01 as an unblended eclipsing 
binary 

With a 3(7 upper limit on the mass of KOI-12.01 of 
about 10Mj up , we can exclude that the photometric 
signature was caused by the transit of an unblended 
low-mass star or brown-dwarf. 

2. KOI-12 blending the light of an eclipsing binary 

The transit light curve could have been produced 
by the occultations of an eclipsing binary in the 
background/foreground, diluted by the light of the 
fast-rotator KOI-12. This scenario can however be 
excluded with tomography, as the time-series spectra 
show that the transit signature is located within the 
rotationally-broadened profile of KOI-12, and corre¬ 
spond to the transit of KOI-12.01 in front of KOI-12. 

3. KOI-12 and KOI-12.01 as an eclipsing binary, blended 
with the light of a tertiary star 

Using PASTIS, we found that the only triple system 
consistent with the Kepler photometry would need KOI- 
12 . 01 to be a small star with Rkoi-12.oi/-Rkoi-12 ~0.1 
and a mass of about 0.14 Mq. The high mass of KOI- 
12.01 would induce a large reflex motion on KOI-12, 
which would result in a large Doppler shift of its line 
profile. In this triple system, the line profile of KOI- 
12 is diluted by the light of the tertiary star. There 
are only two possibilities for the light of this star to be 
undetectable in the observed spectra. In the first case, 
the tertiary star and KOI-12 would need to have simi¬ 
lar projected rotational velocities, radial velocities, and 
brightness. This case is not only highly unlikely, but the 
constraints on the photometry light curve from PASTIS 
favor the second case, with a faint tertiary star. PASTIS 
indicates that KOI-12 must be about five times brighter 
than the tertiary star, which means that the CCF would 
be dominated by the large rotationnally-broadened line 
profile of KOI-12. In this case, any faint dilution from 
the putative tertiary would not be able to mask large 
radial velocity variations of KOI-12 inferred from the 
CCF. Since no significant keplerian motion was detected 
in the data, we can exclude the triple system scenario. 
We note that a binary system in which KOI-12.01 would 
still be a planet can be excluded for similar reasons: 
either KOI-12 and the tertiary star would have to be 
uncannily similar, or the tertiary star would have to be 
very faint and would not impact the measurements. 

We thus conclude to the planetary nature of the Kepler 
candidate KOI-12.01, herafter designated KOI-12b. 

8. Discussion and conclusion 

We assess the existence of a giant exoplanet around 
the Kepler target star KOI-12, with an upper mass 
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Fig. 9. Correlation diagrams for the probability distributions of the tomographic model parameters, without priors from 
photometry. The plot layout is the same as in Fig. [7] 

limit of 10A:fj U p. The inflated radius of K0I-12b (R p = 

1.43±0.13i?j up ), unexpected for this moderately irradi¬ 
ated warm Jupiter (a~0.14au, T eq ~1110K), makes it 
the largest exoplanet known at orbital distance greater 
than 0.1 au. This raises questions about the origin of this 
planetary system, which can be studied through the mea¬ 
surement of its obliquity. We used line-profile tomography 
of SOPHIE time-series spectra to identify the prograde, 

low-obliquity orbit of K0I-12b with A = 12.5-2'2 . Doppler 
tomography allowed a higher precision on the obliquity 
to be reached by comparison with the analysis of the RM 
anomaly detected in the RV measurements derived from 
the SOPHIE spectra, although the two methods yielded 
consistent results. This technique can also be used to 
detect flux anomalies in the spectra, that can bias the RV 
measurements. 

Time-series CCFs used in tomography are generally 
computed using over-sampling with velocity/wavelength 


bins much smaller than the instrumental resolution. This 
results in correlated noise preventing the use of y 2 statis¬ 
tics to adjust the model parameters . In previous studie s 
(e.g. Collier Cameron et al. 2010a Miller et al. 2010), 
uncorrelated Gaussian noise was recovered by resampling 
spectra at the instrumental resolution. Here we devised an 
empirical method, based on the analysis of the tomography 
residuals, to estimate statistically independent error bars 
on the unbinned spectra. Because of the size of KOI-12b 
and the very high rotational broadening of its host star, 
the velocity width of the missing starlight was large enough 
to be detected in the spectra binned at the instrumental 
resolution. In the future, our method should allow for 
smaller planets, or planets transiting stars with very low 
rotational broadening, to be more easily detected using 
tomography. 

As can be seen in Fig. [lOj most obliquity measurements 
have been made for massive exoplanets at close orbital 
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distance and until now only five planets at more than 0.1 au 
(or P ^11 days) had their alignment known. HD80606b 
and HD 17156b are two giant planets on a highly eccentric 
orbit. The former orbits alone its host-star HD 80606, 
in a binary system with HD 80607, and shows a large 


misalignement (iHebrard et al. 2010 

); the latter displays a 

low obliquity (Bartiieri et al.[2009 

Marita et al. 

2009 

and 

possibly has a second planet (Short et al. 

2008 

. The last 


and two smaller plan ets on a coplaner, well-aligned orbit 
(Sanchis-Ojeda et al. 2012). By comparison, KOI-12b is 
an isolated giant planet on a slightly misaligned orbit. It 
is hazardous to search for trends in this limited sample, 
but for now obliquities of planetary systems beyond 0.1 au 
seem lower and less varied than for closer-in planets (see 
Fig. 10). 

With an effective temperature T o ff=6820±120 K and 
a projected rotational velocity 'i>sim*=60.0±0.9kms _1 , 
KOI-12 is one of the hottest star known to host an exo¬ 
planet and the fourth fas test ro tato r h ost afte r WAS P-3 3 
(usini*=86.1±0.4 km s _ 1 , Collier Cameron et al. 2010b I, 
j 73±0.5kms -i , Bieryla et al. 2015 ), and 
70kms 1 


Santerne 


Kelt-7b (i>sim* 

the binary s ystem KOI-13 (vsmi 

|et al.|[201 2b). With an obliquity significantly higher than 
0°, KOI-12b can be considered as slightly misaligned. Yet 
with |A| < 30° the system is no t misaligned in the sense 
defined by Winn et al. (2010a), and does not follow the 
apparent trend that misal igned planets are fou nd around 
hot stars (T e ff >6250 K). |Winn et al. (f2010a|) suggested 
that the thin convective envelop of a hot star limits tidal 
dissipation, preventing the realignement of misaligned 
systems. It is possible that the KOI-12 system acquired 
early-on a moderate obliquity, that it kept to this day 
because of the limited tidal damping. In addition, the 
large orbital distance of KOI-12b and the large mass 
of its host star (1.452±0.093 may have increa sed 

the tidal-dissipation timescale (Albrecht et al. 2012) in 
this relatively young system of 1.5±0.5Gyr. The KOI-12 
system would make an interesting target to study the 
link between obliquit y and (elliptical) tidal instability 
(Cebron et al. 2011), and with a moderately bright 
host-star (Kepler magnitude 11.4) KOI-12b will be a pre¬ 
cious target for future exoplanet atmosphere investigations. 
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